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Gravity beyond linear perturbations in the braneworld 

Hideaki Kudoh*) 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

In the present paper, we briefly review recent studies of second-order gravitational per- 
turbations in braneworld models. After we consider the possibility of pathological behav- 
ior of gravity at higher orders of perturbation, second-order perturbations are discussed in 
Randall- Sundrum braneworld models. Because the mass spectra in these braneworld models 
are different, we analyze them using different approaches. In the respective models, it is 
discussed that 4D Einstein gravity is approximately recovered at level of second-order per- 
. turbation, although there are some exceptional cases in the two-brane model. 
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§1. Introduction 

(N 

The idea of the braneworld has attracted considerable attention from many 
points of view. 1 )" 3 ) In particular, the models introduced by Randall and Sundrum 
Q \ (RS) 2 )' 3 ) have been studied by many people because the so-called warped compact- 

' ification of RS models make the models intriguing. 

A significant point in discussing the viability of these models is whether 4D 
Einstein gravity is recovered on the brane from 5D Einstein gravity in the bulk. In 
the RS models, the gravity that is induced on a flat brane by matter fields confined 
Q" 1 there has the form 4 ) 

6&. ^ Cr> = 8ttU 1 + 



where 8ttG and 87rCr5(:= k) are the 4D and 5D effective gravitational coupling con- 
^ ■ stants, respectively, and tt^ is a tensor quadratic in the 4D energy-momentum tensor 

T^ u . The projected Weyl tensor E^ u is determined by solving the 5D equations. Ap- 
parently, the last two terms are contributions from the existence of the 5D bulk 
spacetime. However, the evaluation of in the above formulation is not easy in 
general, and explicit and typical contributions in general situations have not been 
known, except in the case of linear perturbations. 5 ) One fundamental but important 
fact owing to the warped compactification is that in the RS one-brane model, lin- 
earized 4D Einstein gravity is derived on the brane from 5D Einstein gravity when 
we consider non-relativistic matter on the brane. 6 ) " 9 ) Furthermore, the conventional 
FRW universe is realized on the brane at low energies. These facts imply that the 
last two terms in Eq. (1-1) can be omitted in the restricted configurations. 

The recovery of linearized 4D Einstein gravity in the RS one-brane model is 
expected, and it can be understood by considering the equation for linearized grav- 
itational perturbations, which has a volcano potential due to the tension of the 
brane. 3 ) The volcano potential admits a normalizable bound state, which is called a 
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"zero mode" because the mass eigenvalue of the 4D d'Alembertian is zero, and also 
continuum modes that correspond to massive Kaluza-Klein modes (KK). However, 
the contribution of massive modes to the induced gravity is small compared with 
that of the zero mode, owing to the potential barrier, and therefore the zero mode 
gives the long-ranged force of conventional Newton's law. 

For gravity beyond linear perturbations, the above argument does not apply, 
and for this reason the recovery of 4D linearized Einstein gravity is not expected 
to be realized in a simple treatment. In the one-brane model, if we naively write 
down the 4D effective action, integrating over the extra dimension with the aid of 
the decomposition with respect to the mass eigenvalues, it seems that the KK mode 
coupling is ill-defined and diverge. 10 ) However, this pathological behavior might be 
an artifact of the fact that this argument is based on a mode-by-mode analysis using 
the decomposition of the mass eigenmodes. In fact, the mode-by-mode analysis does 
not give a good description of the behavior of massive modes. 10 )' n ) In the two-brane 
model, other pathological behavior seems to appear on the negative tension brane; 
some hierarchically enhanced terms appear at second order in the perturbation, and 
then seem to prohibit the recovery of the conventional 4D gravity in the nonlinear 
regime. 

Nevertheless, there is the indication that the induced gravity is described well 
by 4D Einstein gravity, even in the highly nonlinear and non-perturbative regimes, 
12 ) -16 ) and this supports the viability of braneworld models. 

In this paper, we review the perturbative approach applied to understanding 
gravity beyond linear-order perturbations, i.e., second-order perturbations. 13 )~ 16 ) In 
the analyses of second-order perturbations in the RS models, static and axisymmetric 
configurations, in which the metric on the branes is spherically symmetric, have 
been assumed for simplicity. In §2 we review the second-order perturbations in the 
RS one-brane model without a bulk scalar field. In §3 we study the second-order 
perturbations in the RS two-brane model where a stabilization mechanism is taken 
into account. Section 4 is devoted to a summary. 

§2. one-brane model 

2.1. Leading -order perturbations 

In this section we review second-order perturbations in the RS one-brane model. 
14 ) The system of this model is simpler than that of the two-brane model, and 
therefore it is more tractable in the study of gravitational perturbations. 

The most general axial-symmetric and static metric has the following form: 

ds 2 = a 2 (y) \- e A ^y)dt 2 + e B ^dr 2 + e c <™VW] + dy 2 . (2-1) 

Here a(y) is a warp factor and it is taken to be a = e~^^, since we consider 
gravitational perturbations around the RS background solution of AdSs, in which 
the vacuum brane of positive tension is located at y = y + = 0. Non-relativistic 
matter fields confined on the brane are assumed to have an energy-momentum tensor 
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of the perfect fluid form 



diag{-p+(r), P + (r), P+(r), P+(r)}. 



(2-2) 



To obtain perturbative equations for the metric functions, we expand A, B, and 
C to second order as A = Ylj=i 2 etc. The 5D Einstein equations with the 
transverse-traceless (TT) conditions 

1 



A W = _ J_ dr (r 3 BM), C« = -d r (r 2 B^) 



2r 

give the master equation for the second order perturbations 



1 



d v a A d v + A 



A (J) 



f a 2 Q yy dy + \d r (r 2 A^B^ 



(2-3) 



(2-4) 



2 + (B :y ) 2 + 2(C, y ) 2 ] . Moreover, we have 
= and e^ 2 ) = 1, to represent the first- and 



where A := Y^=i d i and Qyy = \ [(Ay 
introduced the symbol e^ J \ defined by 
second-order equations in a single expression. The source terms on the right-hand 
side are absent at linear order, because they are quadratic in linear-order quantities. 
The other metric functions at second order, B^ and C^ 2 \ are explicitly related to 
A^ by constraint equations that are obtained from the 5D Einstein equations. 

To solve the master equation, we need to specify the boundary conditions. In 
the original coordinates (2-1), the surface of constant y has been chosen so that the 
metric functions vanish in the limit y — > 00, to coincide asymptotically with AdSs. 
After fixing the surface of constant y at infinity, the coordinates are extended to 
the region near the brane. Therefore, the location of the brane does not generally 
coincide with the y = y + surface. 6 ) In such coordinates, the junction condition is 
not trivial. Thus it is convenient to introduce Gaussian normal coordinates x a , 
in which the location of the brane remains at y = 0, but the metric form is kept 
diagonal. The boundary conditions in the original coordinates are obtained through 
gauge transformations defined by x a = x a + £ a . Israel's junction condition gives the 
boundary conditions as 

2 



v \v=o 



P + + K - 1 (5Af)\y = Q 



(2-5) 



where 5A^ is defined as the difference between the metric in the original coordinates 
(2-1) and that in the Gaussian normal coordinates x a : 8A^> := A^ — A^ 2 \ 

The formal solution to the master equation (2-4) with the boundary conditions 
is obtained by means of Green's function. We have 

[2k 

Y ^ 

'{r 2 A, r B)^ 2 /■■"' 



a 2 A( J \r,y) = 2j dx' 3 G A \ y , =0 {^p^ + {nP + 5A% 
-2e^ J d 3 x'dy f G A 



and Green's function in the static case is given by 



G A (x,y;x',y') 



-I 



d 3 ke ik( X - X ') 



Na(y) 2 a(y 
k 2 + e 2 



l\2 



+ 



+ / dm 



a Qyy 



dy" 



u m {y)u m (y') 
m 2 + k 2 



(2-6) 



,(2-7) 
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where the u m (y) are mode functions, and m is the 4D mass eigenvalue. The normal- 
ization factor N is defined by 

1 := 2 a 2 dy, (2-8) 

where y_ is set to infinity in the present model. The orthogonality of the mode 
functions is expressed by 

r°° dy / 
dyu m = (m/0), 2/ —~u m u m i=5{m-m'). (2-9) 

Jo a 

The mode decomposition and its meaning are obvious from Green's function. The 
first term in Eq. (2-7) is the contribution from the zero mode (m = 0), and the 
second term corresponds to the propagator due to KK states that have non-zero 
mass eigenvalues (m > 0). 

The leading-order in linear perturbations is obtained from the zero-mode trun- 
cation. Substituting only the zero-mode part into Eq. (2-6) and transforming it into 
Gaussian normal coordinates, we obtain the leading-order in linear perturbations 
induced on the brane: 

= 5™ = « -2#+. (aty = y+) (2-10) 

Here we have represented the result by using the Newton potential, which is defined 
by 

A$ + (r) := 4vrG pj } (r) , (2-11) 

where the 4D Newton's constant is defined by 8ttG := kN. Hereafter, the gauge 
freedom of the radial coordinate is chosen to be the isotropic gauge. Now, the 
evaluation of leading-order terms in second-order perturbations is straightforward. 
Substituting the zero-mode part of Green's function and the first-order quantities, 
we finally obtain, at the leading order, 

AA {2 \r, 0) » 8itG (jf® + 3P} 2) - 2<2>+Z\<2>+) , 
AB^(r,0) = AC (2) (r,0) w SirGp^ + 4<P + A<P + - (<2>+, r ) 2 . (2-12) 
These results agree with those for 4D Einstein gravity. 
2.2. Suppression of the KK mode propagation 

The appearance of conventional 4D gravity can be easily understood using the 
zero-mode truncation of Green's function. Thus our main task is to demonstrate the 
suppression of contributions due to KK modes. To do this, we need to evaluate a 
convolution of Green's functions. Because G\ is composed of a zero mode part Go 
and a KK mode part Gk, the contribution from these terms can be decomposed into 
several pieces, depending on which combination of the three propagators is used. 

We note that the mode couplings that contain two zero modes, such as Gk x 
Gq x Go, vanish because of the orthogonality of mode functions. Thus our concern 
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is the mode couplings that include at least two KK modes. Since the mode- by- mode 
analysis yields pathological behavior even at the level of linear perturbations, it is 
necessary to sum up all the mass eigenvalues to handle the KK mode interactions. 
10 ) Carefully performing such analysis, 14 ) it is found that the leading order of the 
corrections is 




(2-13) 



where r± is a typical length scale of perturbations. Therefore, the corrections due 
to the KK modes are suppressed on scales sufficiently large compared to the 5D 
curvature scale I. Note that this suppression essentially arises from the following 
inequality involving Green's function: 

< -G A ( Xl ,y;x' 2 ,0) < _^-(£ + fl 12 ), 

f°° dy l 

/ —G K (x 1 ,y; x 2 , 0)G K (x 3 , y; x 4 , 0) < — T^~^' 

JO a 2 4(4^)^^12^34(^12 + «34j(« + Rl2 + R34) 

where Rab = \&A — xb\- 



§3. Two-brane model 



3.1. Master equations 

In this section, we consider second-order perturbations in the RS two-brane 
model. 15 )' 16 ) The second brane, which has negative tension, is located at y = y_ (> 
y+). Hence, there is a new characteristic length scale of this system, the distance 
between the two branes, which we refer to as the "radius". Because this radius is 
related to the hierarchy, 2 ) it must be stabilized at an appropriate value. As a radius 
stabilization mechanism, we introduce a bulk scalar field, following Goldberger and 
Wise. 17 ) The Lagrangian for the scalar field is 

The scalar field is expanded up to second-order as <p = 4>o(y) + (p^\r, y) + ^ 2 \r, y), 
where 4>o is the background scalar field configuration. The gravitational perturba- 
tions, however, cause the radius to fluctuate, and we cannot assume the metric (2-1) 
in the bulk. Instead, it is useful to assume the "Newton gauge," 

ds 2 = e 2Y dy 2 + a 2 [-e A ~^dt 2 + e B ~^dr 2 + e c ^r 2 dQ 2 ] . (3-2) 

Here A, B, and C correspond to the TT part (2-3), and tp to the trace part. The 
energy-momentum tensors are given in the perfect fluid form as in Eq. (2-2): 

TJ = a± 4 diagj-pi, P±, P±, P±}. (3-3) 

The warp factor in the definition (3-3) is incorporated for the following reason. In 
the present analysis, we adopt the normalization in which all physical quantities 
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are always measured on the positive tension brane at y = y + . Because a length 
scale is warped by the warp factor, physical quantities on the negative tension brane 
at y = y_, such as p_ and P_, are measured at y = y + as cC 4 p_ and cC 4 P_. 
Another normalization, which uses the proper length to measure physical quantities, 
is explained in §3.4. Each normalization has its own virtues. 

The warp factor experiences a back reaction from the bulk scalar field, and in 
general we cannot assume that the warp factor has a pure AdS form. It is determined 
by solving the background equation, 



*%) s [ - 

As compared with the one-brane model (2-1), we have two additional scalar 
functions in the present case, Y and ip. They are not independent, being related to 
each other by the 5D Einstein equations. Moreover, the perturbation of the scalar 
field (p is also related to Y and ip. The constraint equations are given as 



4> {J) (r,y) 
<P {J) (r, V) 



-J—d y (a 2 Y^) + -JL € (.» 
2K(f>oa 2 2k<Pq 



S v + d y A-% 



(3-5) 



Here and Su> are second-order source terms, which are constructed from the first- 
order quantities. 15 ** Here and hereafter we use second-order source terms, such as 5* 
and §*, without giving their definitions. All of them are defined in Ref. 15). 

The equations for and Y^ are obtained from the 5D Einstein equations. 
The equation for A is basically the same as Eq. (2-4), except for the source terms 
on the right-hand side. The master equation for scalar perturbations is given by 



,2 12 



,2^2 



a 2 2 ;2 

o v a —a A - 



(3-6) 



Let us consider the junction conditions for the scalar field. Integrating the equation 
of motion for the scalar field across the brane, we obtain the junction condition in 
the Newton gauge as 



(J) _ 



where we have defined 



(-0 

> il) 



(3-7) 



(3-8) 



Through this boundary condition, the potentials V(±) on the branes affect the per- 
turbations. Combining the junction condition (3-7) and the constraint (3-5), we 
obtain the boundary condition for (3-6), and in the end we can construct the formal 
solution using Green's function. 
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3.2. Derivative expansion 

Formal solutions of the TT part are given by using Green's function (2-7), where 
the continuous mass spectrum in the one-brane model is replaced with the discrete 
one due to the S 1 /Z2 compactification. Because of this difference, we can use another 
method to evaluate contributions from the discrete massive modes. The long-ranged 
contribution of the TT part, Aq, is evaluated using the zero-mode truncation, as 
before. The remaining part, As, which comes from massive KK modes, is evaluated 
using a derivative expansion method. We expand the perturbation variables in terms 
of the expansion parameter (iTr*) -1 , assuming that the typical length scale r+ of 
perturbations is much longer than the 5D curvature scale H —1 . It is important 
to stress that this derivative expansion method is valid only when the mass of the 
first excited mode is sufficiently large. In the one-brane limit (y_ — > 00), the mass 
spectrum becomes continuous and the mass of the first excited mode approaches 
zero (m — >■ 0). Therefore the derivative expansion method is no longer valid in this 
limit. 

For a scalar-type perturbation, there is no zero-mode, owing to the stabiliza- 
tion mechanism. 7 ) To study the contributions from massive modes, we expand the 
perturbation variables in a derivative expansion. The scalar-type perturbations are 
divided into three parts: a pseudo-long-range part Y pse , a short-range part Y$, and 
a part involving interaction terms Y\ . The first one consists of the formal solution of 
the lowest order in the derivative expansion including long-range metric perturba- 
tions. The short-range part consists of the formal solution at the next order, which 
is obtained through iteration of the derivative expansion. The part involving inter- 
action terms, Y\, is defined to consist of the terms that include a potential on the 
brane, represented by A± [see Eq. (3-7)]. They are explicitly given as 

Y^=Y^+Y^+Y^\ 
Mgir, y ± ) = E ^ =F f*lTi J) ~ ^ fa + 2H^Sf) 

<T=± 



-2iVe (J) 



-cr=± Jy+ 2k4>^ 

/ rW t(J)\ 

^(-,«)-^(E^ + | f ), 

4v,, ±) ^r-§£W, ( 3 .9) 

Jy+ Po 

where 



Li J \r) = H(y a ) 



2n4 + a H 2 o ^ un 2^ + 2 K a^3J 



v=v<y 



and the functions u± and v± are given by u± := 1 — 2Hv± and v± := a 2 P a 2 dy'. 

These formal solutions for the TT part and the scalar-type perturbations enable 
us to evaluate the perturbations induced on each brane. 
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To obtain the induced gravity, we must transform the perturbation quantities 
calculated in the Newton gauge into those in Gaussian normal coordinates of the 
form 

ds 2 = d f + fl 2 ( y) (_ e A d p + & B df 2 + e C f 2 dn 2^ (3 . 1Q) 

We need two sets of Gaussian normal coordinates, that in which the positive tension 
brane is located at y = y + and that in which the negative tension brane is located at 
y = y_. After gauge transformations for each coordinate set, we obtain the metric 
perturbations induced on the branes, 



AA { J \r, y ± ) = torG(j%> + 3P± J >) ~ A &A± tV + C ± A ± , r + H{? ± ) 



?S/\2 



+2N^a4(SI-S%) 



+2N / dy 

'y+ 



3u± 



a (Sa± - S,p±) + a v±AS lf ± H r-Z\S y± 



AY, 



(J) 



A ' 



AB( J \r,y ± ) 



AAW k h 
2 T ~Y 



(3-11) 



where the radial gauge of the spatial component at Jth order is the isotropic gauge, 
as before. 

3.3. Evaluation of perturbations 

Now we can explicitly evaluate the gravity induced on the branes. To simplify 
the analysis, we assume that matter fields exist on one of the two branes, which 
is regarded as the visible brane. With this simplification, the sum of the Newton 
potentials #± (:= AttGA^ 1 p±*) which generally appears in the two-brane system as 
we see in Eq. (3-9) is replaced as Yl a =± ^ ~^ ^± a ^ V = V±- 

First, we discuss linear perturbations. Because the evaluation of the spatial 
component is almost the same as that of the temporal component, we henceforth 
concentrate on the temporal component. Assuming weak back-reaction a ss e - ^^, 
the results are summarized as 

4^(r, y + ) = 2d> + + (rlA4> + ) \o (/?+) + O ( 7+ ) - d+ [(^2)' V + «-]}, 

A { l\r, y_) = 2£>_ + O {rlA$ + ) jo (/?_) + O (7-) - P- [aia+ + a_] | , (3-12) 

where a±, f3± and 7± are defined as 
2N*\ ± 2N 3 \± 



a± 



3H± 
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P± ■= -4—0 



e 2 

^, (y = y+) , 

0.1mm\2/l(r 16 \ 4 / £ ■ > ( ' rL,) 



r 

7± == (4^H) _1 = P±(rhs£T 2 



Here, a± is the correction from the interaction terms, the massive mode contributions 
of the TT part are expressed by f3±, and the short-range contributions of the scalar- 
type perturbations are represented by j± , in which we have used the effective radion 
mass fhs- When the radius stabilization mechanism proposed by Goldberger and 
Wise is effective, the mass fhs becomes 0(£ _1 ). 17 - ) 

Let us first consider the corrections induced on the negative tension brane (y = 
y_). Because we are concerned with the hierarchy resolution in this case, 2 ) we have 
set the AdS curvature length t and the hierarchy a + /a_ to the Planck length ip\ 
and 10 16 , respectively. Then the expression for /3_ implies the dominant contribution 
of the massive mode for length scales < 0.1mm. However the KK mode does not 
contribute to the force outside the matter distribution but, rather, to the matter 
energy density, since the gravitational potential appears only in the form A<P±, 
which is proportional to p±. Hence the change of the metric perturbation due to the 
short-range part becomes significant on the negative tension brane only when p^} > 
0(TeV 4 ). This is also true for the correction from 7_ as long as the effective radion 
mass is 0(£ _1 ). For a±, the evaluation depends on the details of the stabilization 
model. This correction becomes important compared to that from the short-range 
part when a 4 _a+ or a_ exceeds 0(l/fn s £ 2 ). We note that when a 4 _a + 3> 1, the mass 
of the stabilized radion becomes small, and then the effect appears also through 7_, 
although this case seems to be exceptional. 

On the positive tension brane, the corrections from (3 + and 7+ are similar to those 
in the one-brane model, and hence are suppressed. As we pointed out, the difference 
is that the corrections in this case are proportional to the matter energy density. 
The factor of a + is suppressed only by H+ / H\ , which is small but not hierarchically 
suppressed. As long as A+ has a natural order of magnitude smaller than £, a + 
is at most 0(-ff 2 _/ij + ). Then, the correction remains less than 0(/3 + )0(r 2 A^ + ). 
However, when A + is much larger than £, the correction becomes larger than that 
in the TT part by a factor of \+/£- Although these choices of parameters are not 
natural, the possibility of an enhanced correction might be interesting. 

Let us discuss second-order perturbations. After straightforward but tedious 
calculations, it is shown on each brane that the leading order in the derivative ex- 
pansion is identical to 4D Einstein gravity even in second-order perturbations. Thus 
to complete the recovery, it is necessary to investigate the contributions from the 
massive modes and the interaction terms, which appear at the next order in the 
derivative expansion. However, we immediately realize that at this order there are 
some terms on the negative tension brane that are hierarchically enhanced by a factor 
of O^fi^/a 2 ^) in comparison to the usual post-Newtonian terms. It therefore seems 
that the recovery of the conventional 4D gravity is spoiled at this order. This point 
is studied in Refs. 15) and 16). As discussed there, however, these dangerous terms 
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are completely canceled out. The remaining contributions are actually suppressed 
apart from some exceptional cases which come from a±, as mentioned above. 
The results are summarized as 

AA^{r,y+) = 8ttG(p^ + 3Pj 2) ) - 4$ + A$ + 

+0[ *A\iom + o(- 



r l J \ \aL J \ a 2 _ 



i( 2 )/V o, \ — a^n/J 2 ) _i_ qp( 2 ) 



+0(P + ) 0(a + ) + 0{ — 



AA K *>(r,y-) = 8irG(p K ' } + 3P y J>) - 4£>_Z\<2>_ 

+0 {^j jo (/?_) + O ( 7 _) + O (J3.) [0(aia + ) + 0(a_)] j. 

(3-14) 

In the case that matter fields are confined to the positive tension brane, the cor- 
rections from KK modes are multiplied by a factor of 0((3+/a 2 _) compared to the 
usual post-Newtonian terms. Nevertheless, the appearance of an enhancement by a 
factor of l/a 2 _ is very likely an artifact of the derivative expansion method, although 
the corrections are still suppressed. Since the condition that the typical length scale 
of the spatial gradient is larger than that of the change in the fifth direction be- 
comes (^ 2 / a - r *) = (fl+/ a ~) ^ 1 near the negative tension brane, P+fa^ appears 
as an expansion parameter. Although the correction seemingly becomes large in the 
y_ — > oo limit, this is due to the limitation of the present approximation. 

3.4. Convention for the length scale 

To end this section, we comment on the convention used for the length scale 
that we have adopted in the previous sections. In our study of the two-brane model, 
we have used the normalization scheme in which the length scale is always measured 
on the positive tension brane at y = y+. Thus, when we regard the negative tension 
brane at y = y_ as the visible brane and investigate physical quantities on this brane, 
we must rescale them to obtain their values on the hidden brane at y = y + . This is 
the reason that we have used the definition of the energy-momentum tensor (3-3). It 
is sometime convenient to measure a length scale r* with respect to its proper length 
scale a(y)r* at the location y. This scheme is adopted in Ref. 7). 

To use the proper length normalization scheme in our analysis, we must change 
some notation and definitions. First, the energy- momentum tensors are defined as 

T ± j? = diag{-p ± ,P ± , P±,P±}, (3-15) 

where and are physical quantities measured on the positive and the negative 
tension branes, respectively. The induced 4D Newton's constant and the Newton 
potential are defined on each brane as 

8vrG± := nNa 2 ± , (3-16) 



A$±{r) := 4ttG±p ( 1 ) . (3-17) 
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The Laplacian operator must accompany the warp factor, i.e., a± 2 A. In spite of 

these changes, the results of the linear perturbations are the same, A { q ] = 2<P±. The 
leading terms of the second-order perturbations (3T4) are given in this normalization 

as 



^AM(r,y±) = 87rG ± (p ! £ ) +3Pl 2) )-4<I>j4r)$±- (3-18) 



The modification just consists of the extra warp factor for the Laplacian and the 
4D Newton's constant. Therefore the analyses presented in previous sections can be 
applied here too with only slight changes. 

However, it is necessary to modify the interpretation of gravitational coupling 
scale. Recall that when we set y = y_ (y = y + ), the negative (positive) tension 
brane is regarded of as the visible brane, and we ignore the matter fields on the 
other, hidden brane. When the negative tension brane is the visible brane, the 
hierarchy resolution is involved, and then the fundamental scale is set to the TeV 
scale, as in the original RS two-brane model. On the other hand, when the positive 
tension brane is regarded as the observable brane, we take the fundamental scale as 
the Planck scale, as in the RS one-brane model. Because the gravitational coupling 
scale on each brane is warped, as described by Eq. (3-16), we must vary the value of 
kN according to which brane is visible. To obtain the corresponding gravitational 
coupling scale at y = y + or y = y„, the modification is given as 



kN 



( 0(M P1 2 ) (y = y+), 

(3-19) 

k O(M^) (y = y-), 



where Mp\ is the 4D Planck mass and Mew is a mass scale of order TeV. With this 
modification, we must change the suppression factors. They are given by 



% (y = v + ) 

/o.imm\2/nr 16 \2/ £ew 2 ( 3-2 °) 

, fed {—) U-2Ttev-0 ' {v=y ) 



where i^>\ ~ M p ^ and fgw ~ ^ew- With these modifications, the analysis given in 
the previous sections for the recovery of 4D Einstein gravity is valid. 



§4. Summary 



We have considered second-order gravitational perturbations in the RS one- 
brane model and in the RS two-brane model with the radius stabilization mecha- 
nism. It is shown in each model that second-order perturbations behave well, and 
the results are basically consistent with 4D Einstein gravity, whose temporal compo- 
nent at second order and spatial component at first order have been observationally 
confirmed to an accuracy of about 0.1 %. 18 ^ 
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In the RS one-brane model, deviations appear as a change of Newton's law 
and are suppressed by a factor of 0((£ 2 /r 2 ) m(r*/£)). Thus it is not feasible to 
observationally distinguish these corrections. 

For the RS two-brane model, the analysis becomes slightly complicated since, 
owing to the stabilized radius, there are scalar-type perturbations in addition to 
the transverse-traceless perturbations. As a model of radius stabilization, we have 
assumed a scalar field that has a potential in the bulk and a potential on each brane. 
Three types of corrections are found: massive mode contributions of transverse- 
traceless perturbations, massive mode contributions of scalar-type perturbations, 
and corrections from the interaction terms on the brane. 

When we consider the case in which the matter fields are on the negative ten- 
sion brane, the correction to 4D Einstein gravity appears at a relative order of 
0((a + /a_) (£/r+) 2 ). With the choice of the hierarchy (3T3), the correction to the 
metric in the linear perturbations becomes comparable to the usual Newtonian po- 
tential when < 0.1mm. However, this correction does not give a contribution to 
the force outside the matter distribution. Hence, it seems that this correction dose 
not prevent reproductions of the predictions of 4D Einstein gravity. We have not 
confirmed if this feature remains in second-order perturbations. Nevertheless, the 
corrections are suppressed by the above factor compared to the usual post-Newtonian 
terms, and therefore the effect due to these corrections is almost impossible to detect 
by observing (post-) Newton's law. 

In the case that the matter fields are on the positive tension brane, the correc- 
tion to 4D Einstein gravity given by linear-order perturbations appears at a rela- 
tive order of 0((£/r*) 2 ). The corrections given by second-order perturbations are 
0((a+/a_) 2 (£/rv) 2 ) compared to the usual post-Newtonian terms. Although the 
corrections are still small, the unexpected extra factor of (a + /a_) 2 might suggest 
that deviations from 4D Einstein gravity appear in higher-order perturbations. How- 
ever, this is very likely to be an artifact due to the limitation of our approximation 
scheme. 

To give a complete proof of the recovery of 4D Einstein gravity, further extension 
of the present analysis is necessary, and a clearer understanding of the recovery 
mechanism in the regime beyond linear perturbations is also necessary. 
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